function y = funcKxR1(nx, nz, zeta, CONSTS)

    [KxS1 ~] = funcKxS(nx, nz, zeta, CONSTS);

    Kx = funcKxfull(nx, nz, zeta, CONSTS);
    
    y = Kx - KxS1;

end